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ABSTRACT 

According to the collapsar model, 7-ray bursts are thought to be produced 
in shocks that occur after the relativistic jet has broken free from the stellar 
envelope. If the mass density of the coUimated outflow is less than that of 
the stellar envelope, the jet will then be surrounded by a cocoon of relativistic 
plasma. This material would itself be able to escape along the direction of least 
resistance, which is likely to be the rotation axis of the stellar progenitor, 
and accelerate in approximately the same way as an impulsive fireball. We 
discuss how the properties of the stellar envelope have a decisive effect on the 
appearance of a cocoon propagating through it. The relativistic material that 
accumulated in the cocoon would have enough kinetic energy to substantially 
alter the structure of the relativistic outflow, if not in fact provide much of 
the observed explosive power. Shock waves within this plasma can produce 
7-ray and X-ray transients, in addition to the standard afterglow emission 
that would arise from the deceleration shock of the cocoon fireball. 

Key words: gamma-rays: bursts - stars: supernovae ~ X-rays: sources - 
Hydrodynamics 



1 INTRODUCTION 

CoUimated flows of plasma with velocities close to the speed of light, commonly referred to 
as relativistic jets, have been discovered in a number of astronomical systems. Objects known 
or suspected to produce them include: extragalactic radio sources (Begelman, Blandford & 
Rees 1984); microquasars (Mirabel & Rodriguez 1999); supernovae (Khokhlov et al. 1999); 
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and 7-ray bursts (GRBs). While extragalactic radio sources produce by far the largest and 
most energetic jets in the Universe, GRBs provide perhaps the most extreme example of 
relativistic flow that may be coUimated, exhibiting speeds of F ^ 100 or more (e.g. Waxman, 
Frail & Kulkarni 1998; Wang, Dai & Lu 2000; Panaitescu & Kumar 2002; Soderberg & 
Ramirez- Ruiz 2002). 

Given the twin requirements of enormous energy ^ 10^^ erg and association with star forming 
regions (see Meszaros 2001 for a recent review), the currently favoured models all involve 
massive, collapsing stars and their byproducts, especially black holes. A "collapsar" forms 
when the evolved core of a massive star collapses to a black hole, either by fallback or because 
the iron core fails to produce an outgoing shock (Woosley 1993; MacFadyen & Woosley 
1999). The shocks responsible for producing the 7-rays must arise after the relativistic jet 
has broken free from the stellar progenitor, whose density is reduced along the rotation axis 
due to an early phase of accretion[| (MacFadyen & Woosley 1999; Aloy et al. 2000; Wheeler 
et al. 2000; MacFadyen, Woosley & Heger 2001; Wheeler, Meier & Wilson 2002). While 
the light, relativistic jet (i.e. light compared to the stellar density) makes its way out of the 
progenitor star, its rate of advance is slowed down and most of the energy output during that 
period is deposited into a cocoon or "wastebasket" surrounding it. The jet head propagates 
at mildly relativistic velocity until it emerges from the edge of the He core into the low 
density H envelope at r^, ^ 10^^ cm (MacFadyen & Woosley 1999; Meszaros & Rees 2001; 
Matzner 2002). At this stage (and provided that the density of the H envelope varies steeply 
with radius; see §3), the head of the jet will advance relativistically, and the cocoon plasma 
could escape swiftly from the stellar cavity and accelerate in approximately the same way as 
an impulsive fireball - its energy will be converted via adiabatic expansion into bulk kinetic 
energy. 

Here we describe the evolution and collimation of such cocoon fireballs. We argue that 
an understanding of the structure and time-dependence of the cocoon plasma can come only 
through a knowledge of the properties of the stellar material through which it propagates. 
A cocoon fireball may be stalled while propagating through an extended or a high mass 
envelope, but could expand freely beyond the stellar cavity of a Helium post- Wolf- Rayet 
star. We show that even if only a small fraction of the energy in the jet (~ 1%) is deposited 



^ The majority of stellar progenitors, with the exception of some very compact stars, will not collapse entirely during the 
typical duration of a GRB. A stellar envelope will thus remain to impede the advance of the jet (see Matzner 2002). 
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into the cocoon, it would have enough kinetic energy to substantially alter the structure 
of the expanding outflow, if not in fact provide much of the observed power. The latter 
may be true for an observer that lies off-axis to the jet or at early times when the initial 
jet contribution to the emission is negligible. We examine the physical conditions within 
the cocoon plasma - namely its confinement, coUimation and initial entropy - to determine 
whether or not this delayed fireball will become matter dominated before it becomes optically 
thin. We discuss the possible role of cocoon fireballs in producing 7-ray and X-ray transients 
(both thermal and non-thermal), along with the standard afterglow emission that would 
originate from the deceleration shock. We suggest that detailed observations of this prompt 
burst and afterglow emission may provide a potential tool for diagnosing the size of the 
cocoon cavity and the initial energy to mass ratio rj — (E/Mc^). It also provides a means 
for probing the state of the stellar medium through which both the initial jet and cocoon 
propagate. 



2 COCOONS: "WASTEBASKETS" OF RELATIVISTIC PLASMA 

The properties of the stellar envelope have a decisive effect on the appearance of a jet 
propagating through it. The characteristic stellar progenitor structure is that of an evolved 
massive star, with ^ 2Mq Fe core of radius ~ 10^ cm and a ~ 8Mq He core extending out 
to ^ 10^^ cm (MacFadyen et al. 2001). In some cases, a cool H envelope reaches out to 
^ 10^^ cm, while in others the envelope has been largely lost. The post-collapse (radiation- 
dominated) pressure profile out to the edge of the He core drops roughly as pHe oc p^f oc 
over some two decades in radius (Meszaros & Rccs 2001). Beyond the He core, it drops dras- 
tically as ph oc r~~^ since at these distances the pressure profile is still the precollapse one. 
The precollapse density of the presupernova model A25 of MacFadyen et al. (2001) scales 
roughly as p oc r~^, so that the radiation-dominated pressure is p oc r"^ (see Meszaros & 
Rees 2001). 

Suppose that a collimatcd beam has been established. If all the particles in the beam are 
ultra relativistic, then pj = \pjC^ and the sound speed is c,. c/\/3. At a given time, the 
beam will have evacuated a channel out to some location where it impinges on the stellar 
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envelope^ at a "working surface" which itself advances out at speed Vh (see Fig. |l]a). If the 
power in the jet, Lj, is roughly conserved and stationary, then approximating the channel 
as a cylinder of radius r, we balance momentum fluxes at the "working surface" to obtain 

;7^V ~ Pen^Vh vk.Vh^Cs, (1) 

where v is the speed of the beam, penv, the stellar envelope density, and 9j the jet opening 
angle. If the beam consists of relativistic plasma then v ^ c, and relativistic fluid mechanics 
must be used (the reader is referred for further details to the generalised formalism devel- 
oped by Matzner 2002). During propagation in the iron and He cores, the head of the jet 
will propagate with subrelativistic velocities. The energy supplied by the jet exceeds that 
imparted to the swept-up stellar material by a factor ~ c/Vh ^ 1. The surplus (or waste) 
energy must then not accumulate near the "working surface" but be deposited within a 
cocoon surrounding the jet (Fig. |I|; slightly resembling the cocoons that envelop jets of radio 
sources: see Begelman et al. 1984 and Begelman & Cioffi 1992). 

After the jet emerges into the H envelope, the sudden and drastic density drop at the 
outer edges permits the jet head to accelerate to velocities close to the speed of light (Vh ~ c). 
Thus, if it is a general property that the jet becomes relativistic near the boundary of the He 
core, the outer edge of the H envelope is reached in a crossing time ~ th/^cTI as measured 
by an observer along the line of sight. For an H envelope density varying as p// oc with 
/9 ^ 3 (see §3), the outer edge of the star is reached in a crossing time that may matter little 
when compare to the He-core traversal time0. The fraction of relativistic plasma injected 
into the cocoon after V/i ~ c will be much reduced (Meszaros & Rees 2001). The amount of 
energy that accumulated in the cocoon while the jet was advancing subrelativistic is then^ 

= / L,{t)dt ^ ^ ^ 5 X lO^^Klo n,ii^..50 erg, (2) 

JO V\i ' 

where tne ~ is the He-core traversal time, V^, the average speed of the jet head (which 

is about c/2; Aloy et al. 2000), and we adopt the convention Q = 10^ Q^, using cgs units. 
The cocoon expands in the transverse direction with a velocity given by K ~ [Pc/ PcmY^"^ ■ 
The forces driving the cocoon expansion will be effective so long aspc > Penv (Matzner 2002). 

^ Reconversion into random energy occurs at the end of the channel, which is a natural site for particle acceleration (Colgate 
1974; Chevalier 1982; Meszaros & Waxman 2001; Ramirez-Ruiz, MacFadyen & Lazzati 2002). 

^ Many presupernova stars, on the other hand, have density profiles with /3 ~ 2 (Chevalier 1989). For these stellar progenitors, 
the jet may be unable to punch through the stellar envelope (Matzner 2002). 

* A more rigorous estimate of the cocoon energy content, which accounts for both relativistic and non-relativistic jet propaga- 
tion, can be found in equation 10 of Matzner (2002). 
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The cocoon material could in principle reach pressure equilibrium with the external stellar 
gas, but this generally would not happen before the jet head has reached the outer layers 
of the progenitor star. The length of the cocoon region is similar to that of the jet, but its 
breadth is determined by its transverse velocity. A wide jet would create a near-spherical 
cavity, but a narrow jet would advance much faster than the transverse expansion of the 
cocoon (Matzner 2002), so that the cocoon would be "cigar-shaped" (or maybe more like 
an "hourglass" in the case when the external pressure is a steep function of radius; see Fig. 
Ic). 

At the radii ~ r^.. (probably the radius of the He core) where the head of the jet starts to 
advance relativistically, the volume of the material deposited into the cocoon, Acav, is related 
to both jet and cocoon expansion velocities by Acav ~ (7r/3)r^(K/V/j)^ (Matzner 2002). At 
that point in time, the cocoon plasma would itself be able to break out and accelerate. Unlike 
the jet, this cocoon material does not have a relativistic outward motion, although it has a 
relativistic internal sound speed (i.e. similar energy to mass ratio). At first an asymmetric 
bubble (since pressure balance may never be reached if the external pressure falls off much 
steeper than r~^) will be inflated which can expand most rapidly along the rotation axis 
(Fig. 0c) and may eventually escape the stellar progenitor. But it may never expand freely 
unless it escapes into an exponentially decreasing atmosphere with /3 ^ 5 (.i.e. bare He star; 
see §3). 

2.1 Trapping and collimation 

Insofar as the cocoon material and the lower-entropy stellar envelope can be treated as two 
separate fluids (i.e. diffusion and viscosity can be neglected), it is feasible to estimate the cone 
angle, 6c, of the expanding plasma. The stellar envelope, which contains the outflow along 
the axis, has a sufficiently large optical depth Tenv ~ 10^^ that most of the radiation released 
is trapped, and transported into the cavity by bulk flow rather than diffusing outwards. The 
optical depth across the cavity is enormous, even if t] is so high that the baryon density is 
low, because of the thermal pair density (T ^ 20 keV), and the radiation is then well enough 
trapped to justify a fluid treatment. Collective plasma effects and magnetic fields may also 
reduce the effective mean free path. Unless there is violent entrainment, there would not 
be much mixing of baryons from the envelope into the cocoon, so that during the build-up 
its ratio of energy to baryon-rest-mass will be given approximately by rj. If the magnetic 
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field contributes significant to the total energy density (i.e. MHD jet) and has a preferred 
orientation, then the pressure and magnetosonic velocities are of course anisotropic, but the 
dynamics would be essentially the same as for pure radiation. When the pressure at the outer 
edges of the cocoon cavity has halved the flow becomes transonic and the cross-sectional area 
is minimised. In this way a directed nozzle can be established. The channel cross section is 
proportional to the total power discharge and varies inversely with the pressure at the outer 
edges of the cocoon cavity (Blandford & Rees 1974). Unfortunately neither the jet thrusts 
nor the pressures at the outer edges of the cocoon are known well enough to quantitatively 
predict the dimensions of the nozzle. Beyond the nozzle, however, the external pressure 
drops steeply, and the cocoon material expands freely in the transverse direction. The flow 
spreads out over an angle ~ 7j^. If its free expansion starts just outside the nozzle, where 
the Lorentz factor of the cocoon material 7c is only ~ 2, then it will spread over a wide 
angle and will develop into a roughly semi-spherical blast wave. The cocoon fireball will then 
expand with 7c oc r/rcav- where rcav ~ Oc'>^* is the typical, initial dimension that the fireball 
would have if it started out spherically. 

3 THE STATE OF COCOON FIREBALLS 
3.1 A brief overview of the fireball model 

Before turning to the question of how the kinetic energy of the cocoon plasma is converted 
to radiation, it is worthwhile to summarise now the essential features of the generic fireball 
scenario. 

In the so-called standard fireball model (see Piran 1999 for a recent review), it is con- 
jectured that the fireball wind, expelled by a central source of dimension Tq ~ 10"'' — 10^ cm 
(notice than in the case of the cocoon fireball, the rclativistic plasma is confined to a much 
extended cavity), accelerates at small radii such that its Lorentz factor grows linearly with 
radius until the entire fireball energy is converted into kinetic energy at rjro (Cavallo & Rees 
1978; Goodman 1986; Paczyhski 1986; Shemi & Piran 1990). This energy must be converted 
to radiation in an optically-thin region, as the observed bursts are non-thermal. The radius 
of transparency of the ejecta is 



where £4^^ is the isotropic equivalent energy generated by the central site. The inertia of the 




(3) 
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swept-up external matter decelerates the shell ejecta significantly by the time it reaches the 
deceleration radius (Meszaros & Rees 1997; Chevalier & Li 1999): 



where the external medium particle density is n{r) — Ar~^, with s = for a homogeneous 
medium n{r) — Uism, and s = 2 for a wind ejected by the stellar progenitor at a constant 
speed. Given a certain external baryon density n(r), the initial Lorentz factor 77 then strongly 
determines where both internal and external shocks develop (see e.g. Fig. 3 of Ramirez- Ruiz, 
Merloni & Rees 2001). Changes in rj will modify the dynamics of the shock deceleration and 
the manifestations of the afterglow emission. 

3.2 The propagation of cocoon fireballs 

The above summary describes the qualitative features of the generic fireball model. While 
similar scaling laws are also expected for the evolution of the cocoon material, the physical 
conditions within and around this plasma are different, and so appreciable deviations from 
the "standard" evolution are thus likely to occur. This variety of propagation effects can 
substantially modify the emergent radiation, so that important constraints on the nature 
of the source producing them can be obtained simply by detecting these changes. For this 
reason, we now turn to consider the cocoon properties in more detail. 

Consider a homogeneous fireball of energy Ec, total mass Mc initially confined to a 
cocoon cavity. Clearly, since r > 1, the initial fireball will be an opaque sphere in thermal 
equilibrium, characterised by a single temperature: Tc ~ 400gi£'c sV^v^g ^eV in a spherical 
cocoon cavity, or T^. f=i 100gEll^2^7^J,to keV in a "hourglass" cocoon. ^'=11/4 for Tc > mgC^ 
(photons and pairs), and it drops to 1 when Tc <S rUeC^ (only photons). 

When the radiation energy dominates the evolution (i.e. rj ^ 1), the fluid expands 
under its own pressure such that its Lorentz factor grows linearly with radius, ■jc oc r/rcav 
When the fireball has a size of r^^ = V'^cacv, all the internal energy has been converted into 
bulk kinetic energy and the matter coasts asymptotically with a constant Lorentz factor 77. 
This is perhaps true for a cocoon propagating into the exponentially decreasing atmosphere 
of a carbon-oxygen or helium post- Wolf- Rayet star and into the circumstellar environment 
beyond it, as in curve (a) of Figure 2. When there is a remaining H envelope (with pn oc r~'^) 
beyond the He core, the deceleration radius of the cocoon fireball may well be inside the star. 
The energy required to sweep up an external stellar mass of rUenv is 7c"^envC^, where 7c oc r 




(4) 
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and menv oc p{r)r^ oc r~^~^^. Thus, the cocoon blast wave is undecelerated for /3 ^ 5. By the 
same token, the initial relativistic jet will expand freely provided /? ^ 3. For a moderately 
dense H envelope, the cocoon fireball, which starts being decelerated by the stellar matter 
at r < (see curve (b) of Fig. 2), emerges from the H envelope with a Lorentz factor 
7c < T). For an extended or denser H envelope, corresponding to curve (c) in Figure 2, the 
cocoon of relativistic material would be stalled before emerging. However, even in this case, 
the cocoon may have more energy than the binding energy of the envelope, and could give 
rise to a "hypernova" (Iwamoto et al. 1998; Paczyhski 1998; Wang & Wheeler 1998). An 
accompanying GRB will be present depending on whether or not the jet is also choked. 



3.3 The role of rj 

The baryon load build-up in the cocoon reservoir, Mc ~ Ec/rjc^, influences the fireball 
evolution by increasing the opacity and thus delaying the escape of radiation. The pair 
opacity, Tp, decreases exponentially with decreasing local temperature, and falls to unity 
when Tp w 20 keV. The matter opacity, on the other hand, drops as r~^, and thus the 
escape temperature may drop far below Tp if > Tp = 1 (Shemi & Piran 1990; Piran 1999). 
Two critical values for rj determine the order of these transitions: 



'SalEraT^' 



1/2 



and 

* = ) - 5 X 10^' E^Hr^P:;''. (6) 

yieimpc^ri^^ J 

The effect of the baryons is only negligible when rj > rjp and the evolution is that of a pure 
photon-lepton fireball (r^ = 1 > Tb). If t] is less extreme, there are two qualitative changes 
in the fireball's mode of propagation. First, the matter opacity becomes important when 
rjp > rj > rjh ^ 5 X 10^ -Si,52^cav(9^c,o^"^- The comoving temperature, in this case, decreases far 
below Tp before t — reaches unity, yet the fireball continues to be radiation dominated 
and most of the energy still escapes as radiation. Second, for rj smaller than rjb, the fireball 
becomes matter dominated before it becomes optically thin, and most of the initial energy 
its converted into bulk kinetic energy (this is likely to be the common situation for the initial 
jet fireball; see §3.1). These two modes of propagation will be referred to in the following as 
radiation dominated and matter dominated, respectively. 
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4 MASSIVE PROGENITORS, STALLED COCOONS AND FE LINES 

The properties of the stellar envelope determined in the previous section have an important 
effect on the appearance of a jet propagating through it. While the details remain uncertain, 
preliminary calculations suggest that a relativistic jet can be launched along the progenitor 
rotation axis (MacFadycn & Woosley 1999; Aloy et al. 2000; MacFadycn, Wooslcy & Hcger 
2001). One would expect baryon contamination to be the lowest near this axis, because 
angular momentum flings material away from it and material with low-angular momentum 
falls into the black hole (Fryer & Heger 2000). The jet would be expected to break free of the 
H envelope, and in principle lead to a successful GRB, provided the central engine feeding 
time exceeds the He-core crossing time and that the stellar pressure continues to drops at 
least as fast as pn oc r^"^. In addition to the stellar pressure required to stall the initial jet, 
there is another, smaller pressure which can prevent the cocoon material to expand freely. 
The cocoon fireball will, in general, be easier to confine than the initial jet, for two reasons: 
first, the fraction of matter with a clear line of sight along the rotation axis will be much 
reduced (i.e. fic > ^^funnei ^ ^j), and second, the cocoon blast wave may never expand freely 
for an H envelope density varying as ph oc r~^ with (3 < 5. 

The energy of the relativistic material that accumulated in the cocoon while the jet was 
advancing subrelativistically is nevertheless much larger than the binding energy of the H 
envelope. As soon as the jet penetrates into the low-density envelope beyond r*, the cocoon 
plasma would itself be able to "break-out" and expand through the envelope along the 
direction of least resistance, which is likely to be the rotation axis of the stellar progenitor 
(perhaps further channelled as the jet penetrates further into the envelope). The cocoon 
fireball, in this scenario, may start to decelerate before it becomes matter dominated. It 
starts being relativistic during the acceleration stage, but a self similar phase could then 
begin after enough external material has been collected. In the interim, the Lorentz factor 
will drop faster than that of a cold fireball propagating through a similar density profile 
would (Fig. 2). If the H envelope did not exist, then the cocoon fireball would reach the 
outer edge of this envelope after a few seconds, corresponding to case (a) in Fig. 2, where 
tc^H fa rcav/c+rH/{2cri'^). However, for very extended or slow rotating stars (i.e. denser cores), 
the cocoon material would generally carry less energy and inertia than the stellar envelope; 
it would then take up to a few hours for it to reach the outer surface of the star as it expands 
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subrelativistically with a velocity of the order of Vc ~ c(£'c/MenvC^)^^^ ~ 10^ (Mgnv/M©)"^/^ 
cm (the case of an ultimately choked cocoon). 

As it expands relativistically (but notice that, as argued above, for very extended or slow 
rotating stars the cocoon material may be stalled prior to the acceleration of the fireball 
to relativistic velocities), the cocoon fluid cools with T oc (r/rcav)""*^- The coasting photons, 
whose local energy is T are blue shifted. An observer detects them with a temperature of 
Tohs oc 7c(r)T(r). Seeing that T oc and jc oc r, we find that during the acceleration stage 
^obs ~ Tc, where 

70 r-l{'L]%{^)-y^ keV (spherical) 

(7) 

^4:r~^ii^Ly^oV~g^^ keV (hourglass) 

is the initial black-body temperature of the opaque, cocoon cavity. These photons would 
escape freely after the expanding fireball becomes optically thin, which is likely to occur at 
some distance from the stellar surface and thus Tobs ^ T^. The escape temperature may drop 
far below 7c (r) 7], if condition rj > rji, is not satisfied or if the deceleration of the fireball takes 
place before the radius of transparency (but it is always ^ 40007c(r) K, the recombination 
temperature). A similar diagram to Fig. 2 can be drawn for Tg^s as a function of radius (see 
Fig. 3). 

In the case of a successful break-through of the jet, a strongly decelerated, cocoon fireball 
could result in a potentially interesting and observable phenomenon. Not only would a 
conventional "long" GRB be detectable, followed by a standard afterglow, but also there 
would be, after some seconds or up to a day, a secondary, almost thermal (see Goodman 
1986) brightening , caused by the cocoon photospheric emission containing ^ 10^° erg . The 
observed duration of this emission would be of the order of ?"rc/(c7c) ~ a few hours for 
a cocoon expanding at mildly relativistic velocities (case (b) in Fig. 3). A choked cocoon 
will, however, expand more or less isotropically through the rest of the envelope, causing its 
disruption. This would then appear, after the disrupted envelope becomes optically thin, as 
a type II supernova (or a type Ib/c if there is a significant injection of radioactive material). 
Some evidence for supernova-type emission has been found in: GRB 980326 (Bloom et al. 
1999); GRB 970228 (Reichart 1999); GRB 990712 (Bjdrnsson et al. 2001); GRB 000911 
(Lazzati et al. 2001); GRB 011121 (Bloom et al. 2002; Dado et al. 2002; Garnavich et 
al. 2002). The decelerating external shock of a cocoon fireball may produce a significant 
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absorption edge in the cooled shocked ejecta, so long as cooling is faster than adiabatic 
losses and protons are well coupled to the electrons. Absorption edges can also arise from 
cooler, denser material or filaments in pressure equilibrium with the shocked envelope ejecta 
(Meszaros & Rees 1998). 

It is, of course possible that the initial jet is a magnetically confined configuration, whose 
coUimation properties are unaffected by the distribution of the external matter. In this case, 
the cocoon cavity would still have a dynamically-important magnetic field strength. If the 
field were tangled, continuing rcconnection processes may lead to acceleration of non-thermal 
electrons. It is therefore plausible that a substantial fraction of the energy stored in the 
cocoon cavity could be released, via magnetic dissipation, in a non-thermal ~ UV/X-ray 
continuum with L Ec/tc,H ~ 10^^ erg s~^. A magnetic field of 10^ G could confine clumps 
of gas with densities up to n ^ lO^^cm"^, even at keV temperatures. Such clumps would 
be optically thick and could reprocess a non thermal UV/X-ray continuum arising from 
within the dilute plasma between them (as envisaged by Meszaros & Rees 2001) and (in 
the absence of both entrainment and substantial pair production) may also excite Fe-line 
emission. If the UV/X-ray continuum can maintain a ionization parameter ^ = Lx/{nr'^) of 
the order of 10^ — 10^, a modestly supersolar Fe mass fraction could yield a recombination 
line luminosity comparable to the one observed in GRB 991216 (Piro et al. 2000; Ballantyne 
& Ramirez-Ruiz 2001; Vietri et al. 2001; McLaughhn et al. 2002; Ballantyne et al. 2002; 
Kallman, Meszaros & Rees 2002). 



5 COMPACT PROGENITORS, GRB PRECURSORS, AND AFTERGLOW 
SIGNATURES 

A relativistic cocoon fireball is likely to escape if the star loses its hydrogen envelope before 
collapsing (MacFadyen & Woosley 1999; Aloy et al. 2000; Wheeler et al. 2000; MacFadyen 
et al. 2001; Matzner 2002). This is expected for example in stars with high radiative mass- 
loss (e.g. Ramirez- Ruiz et al. 2001). The fireball cocoon will then emerge from the He core 
into an exponentially decreasing atmosphere and into the rarified circumstellar environment 
beyond it, where it acquires the limiting bulk Lorentz factor 7^ fa 77, unless of course, r] > rji,. 
In this later case, the cocoon plasma continues to be radiation dominated when it becomes 
optical thin. At this stage, the baryons will switch immediately to a coasting phase with 
7p Ri r-r^/rcav (where r^-^ < r^) and most of the energy escapes as photons. As follows from 



12 E. Ramirez-Ruiz, A. Celotti & M. J. Rees 

the previous discussion, an observer will detect them with a temperature of Tobs oc 7cT. 
Thus the observed thermal peak frequency would be in the BATSE [20-600] keV spectral 
window (see equation 7) for compact He envelopes r* ^ 10^^ cm (the outer edge of the He 
core varies with initial mass, roughly as r* ~ 10^^ {Mi/ 35 Mq)'^; Woosley, Langer & Weaver 
1993) . For extended He envelopes (r* ^ 10^^ cm), however, the thermal emission could be 
detectable with instruments like Ginga and the 5e^^oSAX wide field cameras. The time 
delay between the cocoon photospheric emission and the start of the main burst is given by 
At^ fa Tcav/ (2c) +rr^/{2cr]'^) — max [r^-j , Tint]/ (2c77^), where Tint is the radius of internal shocks 
in the jet fireball. The thermal signal emerging from a cocoon fireball with rjp > rj > r)b will 
most likely appear as a transient signal at the beginning of the main burst since | Ai^| ^0.1 
s. In contrast to the main burst, this signal would only last for Tcav/c ~ 0.1 seconds. Such 
thermal precursors may have indeed been observed by Ginga (Murakami et al. 1991). For 
compact stellar cocoons (rcav < 10^° cm), however, this thermal signal could be very short 
lived and may be difficult to disentangle from the internal shock emission. 

The cocoon fireball will strongly modify the usual properties of the standard internal 
shocks and the afterglow emission when rj < rjb (i.e. the cocoon fireball becomes matter 
dominated). The total energy available in the cocoon afterglow is essentially, the kinetic 
energy of the relativistic wind deposited during the He-core traversal time, here estimated 
to be ^ 7r*^ii seconds, minus the fraction dissipated (and converted into prompt 7-rays) 
in internal shocks. The size of the stellar cavity and the initial Lorentz factor 77, along 
with the jet lifetime, strongly determine whether or not the cocoon fireball would carry 
less energy and inertia than the relativistic jet itself. For E^^^^c > E^t^j, the main afterglow 
will be dominated by the deceleration of the cocoon fireball (this scenario is hkely to occur 
when the source hfetime is of the order of the He-core traversal time, so that the energy 
carried by the initial jet is much less than that accumulated in the cocoon). An interesting 
consequence of this scenario is that if the observer lies off-axis to the jet, there could be 
a large fraction of detectable afterglows for which no 7-ray event is detected (commonly 
referred to as"burstless" or "orphan" afterglows; see Meszaros, Rees & Wijers 1998). The 
initial jet starts being decelerated by the external medium at a smaller radius, so that the 
cocoon material always overtakes it and sweeps the jet material. This collision would be 
an important contribution to the observed afterglow at early times (~ < ^^), when 
the afterglow emission produced by the jet fireball dominates. The afterglow shock may 
experience, after starting out in the canonical manner, a "resurgence" similar to the shock 
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refreshment produced by the delayed energy injection of a long-hved central engine (see 
Panaitescu, Meszaros & Rees 1998 for the case of GRB 970508). 

If the jet produced by the accretion maintains its energy for much longer than it takes the 
jet to reach the surface of the H envelope, or is very highly coUimated, the initial jet is likely 
to be more energetic than the cocoon material (i.e. Ej > Ec, and since entrainment may 
dominate in the initial fireball E^-^j = > £'4^,^ = ^jf^)- Figure 4 shows the evolution of 
both jet and cocoon fireballs for this matter dominated case. The adiabatic fireball evolution 
(thick solid line) was computed using a similar numerical method to the one developed by 
Kobayashi, Sari & Piran (1999). The main afterglow will be produced by the slowing down 
of the jet as in the usual case, however, the emission at early stages would be caused by the 
deceleration of the cocoon fireball. Figure 5 shows the contribution of the cocoon afterglow 
radiation at low frequency and at early times when the initial jet contribution to the emission 
is negligible. If the external medium is homogeneous the sub-millimiter afterglow produced 
by the jet should rise slowly at times between a few hours and one day, while for the cocoon 
material the emission should fall steeply after this. Therefore observations made at sub- 
millimiter frequencies with the SCUBA (James Clerk Maxwell Telescope) or with MAMBO 
(IRAM telescope) instruments would be very powerful in determining if the early afterglow 
is dominated by the cocoon fireball emission. 

There could also be additional precursor signatures which are not associated with the 
ejecta blast wave, but with the dynamics of the fireball at the coasting phase. The 7-ray (i.e 
internal shocks) signal emerging from the cocoon fireball would precede by Atr 
few seconds that produced by the initial relativistic wind (provided that E^^^j > E^^^^c ^-nd 
> Tcav/c). The observed variability time scale of this prompt gamma-ray emission is 
related to the typical size of the shocked plasma region containing the photon field: A Tcav 
or A ^ for r > A'j'^. For compact stellar cocoons (i.e. r* < 10^^ cm), the expected delay 
between the cocoon precursor and the main pulse, Atr, would be proportional to the total 
initial jet energy: E^t^j. It will then simply refiect a correlation between the burst strength 
and the time elapsed since the previous emission episode, similar to what is observed in 
GRB lightcurves (Ramirez-Ruiz & Merloni 2001). 
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6 DISCUSSION 



We have discussed, in the context of a collapsar model of gamma-ray bursts (which are 
normally assumed to occur in shocks taking effect after the relativistic jet has broken free 
from the stellar envelope) , the dynamics and evolution of the relativistic plasma surrounding 
a light, relativistic jet. As the jet makes its way out of the stellar envelope, most of its 
energy output during that period goes into a cocoon of relativistic plasma. This material 
subsequently escapes along the direction of least resistance. Provided that the density along 
the H envelope vary monotonically with radius as pn oc r^*^, the properties of the cocoon 
plasma would be similar to those argued in §3. A collimatcd cocoon moving into a region 
with P Ki 5 (i.e. no H envelope) will become overpressured relative to its surroundings, and 
thereafter expands freely. If the relativistic jet carries less energy and inertia than the cocoon 
plasma itself (^^ > %^), it will start to decelerate at a smaller radius than the coUimated 
cocoon fireball, so that the latter would overtake it. In this case, the afterglow would be 
dominated by the emission of the cocoon material, which is likely to be ejected at larger 
angles relative to the observer than those from the jet itself. On the other hand, if the jet 
produced by the accretion maintains its energy for much longer than it takes the jet head 
to reach the surface of the He envelope, the relativistic jet is likely to contain substantially 
more energy than the off-axis cocoon material (since ftj <^ D,c), so that it dominates the 
flux after expanding for a longer time than the initially observed off-axis region. 

Additional effects are expected when the cocoon fireball material becomes optically thin. 
Shock waves within the plasma can contribute with a short-lived (~ few seconds) non- 
thermal 7/X-ray transient for a cocoon propagating inside a compact post- Wolf- Rayet star 
(and so long as 77 < ~ 5 x 10^ El^52f^^^Jg6^ Q^^; see §3). For very extended or slow rotating 
stars, a long lasting (few hours to a day) UV/X-ray (almost) thermal pulse, whose total 
energy may be a few percent of the total burst energy, is likely to appear. If magnetic 
dissipation within this plasma is important, it is also possible that a substantial fraction 
of the energy stored in the cocoon can contribute a non-thermal UV/X-ray afterglow, and 
also excite Fe line emission from the envelope gas. The detection of these prompt multi- 
wavelength signatures would be a test of the collapsar model; and the precise measurement 
of the time delay between emissions may help constrain the dimensions and properties of 
the H envelope, the size of the cocoon cavity, the initial dimensionless entropy of the jet, 77, 
and the radius of the emitting region. The processes discussed here suggest that if GRBs 
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are the outcome of the collapse of massive stars involving a relativistic fireball jet, bursts 
and afterglows may have a more complex spectra and time-structure than those alluded to 
the "standard" model. 
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Figure 1. Schematic diagram illustrating the propagation of a relativistic jet through the stellar envelope. Initially, the jet is 
unable to move the envelope material to a speed comparable to its own and thus is abruptly decelerated (a). Most of the energy 
output during that period is deposited into a cocoon or "wastcbasket" surrounding the jet (b). After the jet head advances 
relativistically, the cocoon plasma would itself be able to escape swiftly from the stellar cavity and accelerate in approximately 
the same way as an impulsive fireball (c). 
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Figure 2. This diagram shows, for three illustrative cases, how the cocoon expansion would be affected by the properties of the 
stellar envelope through which it propagates. The axes (logajrithmic) are 7c versus r, where r is rcav ~ r^dc when the cocoon 
is observed to start its expansion. Three illustrative cases are depicted. In case (a), the stellar matter has a low density (i.e. 
exponentially decreasing atmosphere of a carbon-oxygen or helium post- Wolf- Rayet star), and the cocoon blast wave sweeps 
up all the envelope's matter before it has been decelerated. When the fireball has a size of r^^, = Tjrcavj all the internal energy 
has been converted into bulk kinetic energy and the matter coasts asymptotically with a constant Lorentz factor r/. In case 
(b), with higher stellar density, deceleration occurs at radii < r^^, and the blast wave is still moving through the H envelope 
material during the afterglow. After that, the cocoon escapes into the circumstcUar environment beyond the stellar envelope, 
where it escapes freely with 7c < ri. For a very dense (or extended) H envelope, corresponding to case (c), the cocoon material 
would be unable to break free from the stellar envelope. However, even in this case, the cocoon material has more energy than 
the binding energy of the envelope. 
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Figure 3. Schematic plot of the escape temperature Tobs as a function of radius for the three, qualitatively different cases 
depicted in Figure 2. As it expands, the cocoon fluid cools with T oc (r/rcav)~^- The coasting photons, whose local energy is 
T are blue shifted. An observer would detect them with a temperature of Tobs lc{r)T(r). These photons would escape freely 
after the expanding fireball becomes optically thin, which is likely to occur at some distance from rn- 
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Figure 4. Evolution of a matter dominated, cocoon fireball from its initial formation at rest to its final stage. Both the initial 
jet (solid line) and the cocoon fireball (dashed dotted line) propagate into the exponentially decreasing H atmosphere and 
into the circumstellar environment beyond it, where they will both accelerate while expanding and then coast freely until the 
surrounding matter will eventually influence their coasting expansion. The energy dissipation is due to intera<:tion with the ISM 
via a relativistic forward shock and a Newtonian reverse shock. The parameters for this computation are: r] = 10, rcav ~ 10® 
cm, Lj = 10^" erg s~^, nigm = 1 cm~^, 6j=0.1 and 8c=l-0- Shown are the numerical value of the average Lorentz factor (solid 
and dashed lines) and its analytical estimate (dotted line) 
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Figure 5. Comparison between the initial jet and the cocoon afterglow lightcurves at early times and at observing frequency 
V = 10^^ Hz. The afterglow brightness depends on the relationship between this frequency and those of the injection (i-'m), 
cooling (uk), and absorption (i/q) breaks (e.g. Panaitescu & Kumar 2000). The dotted and dash lines represent the jet and 
cocoon afterglow contributions corresponding to the parameters of the calculations shown in Figure 4. 



